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We formulate a new Hartree-Fock-Bogoliubov method applicable to weakly bound de- 
formed nuclei using the coordinate-space Green's function technique. An emphasis is put 
on treatment of quasiparticle states in the continuum, on which we impose the correct 
boundary condition of the asymptotic out-going wave. We illustrate this method with 
numerical examples. 



1. Introduction 

The Rl-beam facilities in the new generation will increase significantly the number 
of experimentally accessible nuclei, especially in medium and heavy mass regions. 
We may reach nuclei close to the neutron drip- line in the 10 < Z < 20 and TV > 20 
region, where a bunch of deformed neutron-rich nuclei are expectecffi This will 
provide us with a new opportunity to study interplay among the presence of weakly 
bound neutrons, the coupling to nuclear deformation effects, the pairing correlation, 
and the collective excitations. 

A promising theoretical framework to describe this situation may be the self- 
consistent mean-field approaches. ™ More specifically, we consider here the Hartree- 
Fock-Bogoliubov (HFB) method to construct the pair-correlated and deformed 
ground state, and the quasiparticle random phase approximation (QRPA) to de- 
scribe the excitation modes built on the ground state. As commonly recognized, 
one has to describe precisely the nucleon wave function of weakly bound and un- 
bound orbits, which should hav e proper asymptotic behaviours. There exist such 
formulations for spherical nuclei EEEElZl but a new challenge here is that we have 
to do it for deformed nuclei. A method using the Poschel-Teller-Ginocchio basis is 
proposed recently 21 . The quasiparticle motion in deformed Woods-Saxon potential 
is analyzed in detail in the coupled-channel formalism.^ We here take a slightly 
different approach based on the coordinate-space Green's function technique since 
we plan to apply it also to the continuum QRPA^. In the present work, we shall 
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show that the coordinate-space Green's function technique enables us to formulate 
the deformed continuum HFB method in which the nucleon waves satisfy a proper 
boundary condition of the asymptotic out-going wave. 

2. Deformed continuum HFB method using the Green's function 

We first describe the quasi-particle motion in the HFB mean-fields consisting of the 
particle-hole field and the pair field which are both deformed. The axial symmetry 
is assumed. The quasiparticles of the Bogoliubov type have two-component wave 
functions ijj^ 1,2 ^(ra), for which we use the radial coordinate system and the partial 
wave expansion 

L max 
L 

with L = (jlm) and yij{ra) being the spin spherical harmonics. The HFB equation 
is then written as a coupled-channel Schrodinger equations^ for the radial wave 
functions (r)} where the quantum number L represents the "channel". Note 
that the energy spectrum of the quasiparticle consists of discrete and continuum 
parts, which are separated by the energycondition E < |A| and E > |A| (A is the 
Fermi energy) as is in the spherical case!^ 

We can construct the exact Green's function for the quasiparticle motion in 
our deformed HFB problem. It is an extension of the spherical theory of RefP to 
deformed cases, and we accomplished this by employing a general prescriptiorP^ 
of constructing the exact Green's function for a deformed potential scatterer. Here 
the HFB Green's function (a 2 x 2 matrix form combining the normal and abnormal 
functions) is expanded as 

N c 

G(r<T,r'a',E) = £ yL(ra)g LL ,(r, r' , E)y[,(r'a'). (2) 

L,L' 

The coupled-channel radial Green's function gLL'ir, r' , E) is constructed as a linear 
combination of products of "regular solutions" {0^1 (r)} (i.e. those satisfying the 
boundary conditions ^%{r) — > r'5iL"5y at the origin r — » 0) and "out-going wave 
solutions" {^2z/ ( r )} (those connected to the proper asymptotic form 4>^)(r) — ► 
r _1 Hjf (kir)SLL'Sij for r — > oo where H^{kr) is the out-going Hankel function). 

We calculate the density p(r) and the pair density p(r) using the HFB Green's 
function thus constructed. The generalized density matrix 

R(ra,rW) = ( ^'fJ. . f ^') = ± f G(ra, rV, E)dE, 

\p*{ra,r'a') Srr'Saa' - p(ra,r'a') J 2m J c 

(3) 

which is a sum of the wave functions of all the quasiparticle states including the con- 
tinuum states, is calculated using a contour integral of the HFB Green's functional. 
Incorporating this way of calculating densities into the standard iterative algorithm, 
we obtain the HFB ground state after convergence. 
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Fig. 1. The monopolc and quadrupolc parts, po(r) and P2<V)> °f * ne neutron density (left panel). 
The monopole and quadrupole parts, po(r) and p~2(r), of the neutron pair density (right). The 
dashed curves are the results obtained with the box boundary condition. 



3. Numerical analysis 

We shall demonstrate with numerical examples how the deformed continuum HFB 
works. We adopt for simplicity a deformed Woods-Saxon potential as the particle- 
hole field, but we perform the HFB iteration to obtain the sclfconsistent pair field. 
We use the density-dependent delta interaction (DDDI) acting in the singlet pair, 
Vpair = w (l — ?7(/9„(r)/0.08) - 59 ) 5(r — r'), (for neutrons) where vq is fixed to 
reproduce the nn-scattering length a = —18 fm. We consider 38 Mg and assume a 
deformation (3 — 0.3. Using the Runge-Kutta-Nystrom method we solve numerically 
the coupled-channel equation within an interval r = [0, r max ] (r max — 15 fm) with 
a step size Ar = 0.2 fm. At the outer boundary the wave functions are connected 
to the asymptotic forms. The parameter r\ = 0.76 is chosen to produce the neutron 
pairing gap around A ^ 1.5 McV. The cut-off in the quasiparticlc energy is 60 MeV, 
and the maximum fl (j z ) quantum number is fl m ax — 21/2. For comparison, we 
performed also the HFB calculation using the same model but with a box boundary 
condition assuming an infinite wall at r = r max . In the following we show results 
for neutrons. 

Figure 1 shows the radial profile of the monopole and quadrupole parts po{r) and 
P2{r) of the neutron density p(r) = J2\ PaM^ao^), and the corresponding po(r) 
and p2{r) of the neutron pair density p(r). We obtain exponential asymptotics here 
thanks to the proper boundary condition, and it is in contrast to the results obtained 
with the box boundary condition (the dashed curves in Fig.l). po(r) and p2(r) have 
the same exponential slope, indicating that we can define the deformation of the 
equi-density surfaces in the asymptotic region. This kind of deformed exponential 
tail is also seen in the neutron pair density. But the ratio of p~2{r) against po{r) is 
significantly smaller than that of the normal density. This points to that the pair 
density in the tail has smaller deformation than that of the normal density. 

The quasiparticle spectrum above the threshold energy E t h = |A| should be con- 
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Fig. 2. The occupation number density n(E) (left panel) and the pair number density n(E) (right 
panel) for neutrons, plotted with the solid curves. The results obtained with the box boundary 
condition are also plotted with the dashed curves. The inset is a magnification of n(E), and 
we compare it with the result (the dotted curve) obtained with a Woods-Saxon potential whose 
bottom is shifted up by +2 MeV. 



tinuous, and it is indeed the case in our formulation. Figure 2 shows the occupation 
number density n{E) and the pair number density h{E) which are defined by 

n(E) = —lm^2 J drGu(ra,ra,—E—ie), h(E) = — Im^ J drGi2(ra,ra, -E-i 

(4) 

They quantify contributions of the quasiparticle state at energy E to the neutron 
number J p(r)dr = N and to J p{r)dr. Here is shown only for the neutron CI = 1/2 
states. The smoothing parameter is chosen to e = 25 keV (a discrete state would 
have artificial FWHM=50 keV). 

In this numerical example there is no discrete quasiparticle states below |A| (A = 
—889 keV), and all the quasiparticle states are embedded in the continuum E > |A|. 
The quasiparticle states corresponding to deep hole Woods-Saxon orbit appear both 
in n(E) and n(E) as narrow resonances. It is observed also that the non-resonant 
continuum states have a significant contribution to the pair number density h(E). 
This figure shows also that the box-discretized calculation has difficulty to describe 
the non-resonant continuum states. The lowest energy resonance is not described 
well by a single state in the box-discretized calculation, and it is because this 
resonance has a rather large width. This resonance corresponds to the [310]^ orbit 
in the deformed Woods-Saxon potential, which is, in the absence of the pairing, a 
bound state with the single-particle energy e = —798 keV. 

Naturally we expect most dramatic effect of the weak binding on this state. 
The inset shows how the spectrum h(E) changes when the depth of the Woods- 
Saxon potential is artificially shifted (made shallower) by 2 MeV. The Woods-Saxon 
single-particle energy e and the Fermi energy A changes from e, A = —798, —889 
keV to e, A = —40, —88 keV. We see in the inset of Fig. 2 a dramatic increase in 
the width of the lowest-energy resonance, which apparently originates from the 



6, 2009 9:17 WSPC/INSTRUCTION FILE JFS 



Continuum Coupling and Pair Correlation in Weakly Bound Deformed Nuclei 5 

weak binding. Note however that the peak energy of the resonance stays almost 
constant. This implies that the effective pairing gap of this resonant quasiparticle 
state is unchanged, if we estimate the effective pairing gap through the relation 
E qp = y/{e — A) 2 + A 2 . This observation is different from that in RefP^ claiming a 
reduction of the effective pairing gap due to the weak binding effect. The difference 
originates from the fact that we here use the selfconsistent pairing field generated 
from the DDDI, whose force strength becomes large at low densities. The pair field 
extending to far outside the nucleus plays important roles. 

4. Conclusions 

We have formulated the deformed continuum HFB method which is designed for 
weakly bound deformed nuclei. We utilized here the exact construction of the quasi- 
particle Green's function for deformed HFB mean-fields on the basis of the coupled- 
channel representation. The proper boundary condition of the out-going wave is 
imposed on the continuum quasiparticles. We have analyzed numerically effects of 
the continuum coupling and the weak binding on the pair correlation. It is found 
that the quasiparticle states in the non-resonant continuum play significant role to 
generate the pair correlation. It is also suggested that the effective pairing gap of 
weakly bound orbits is not reduced very much, provided that the pairing interaction 
has the surface enhancement. 
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